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ABSTRACT 

If  a  weapon  system  is  fired  at  a  circular  target  and  the 
impact  point  is  distributed  as  a  bivariate  normal  random 
variable,  it  is  not  particularly  difficult  to  determine  hit 
probability  if  the  expected  impact  point  is  at  target  center. 
If,  for  some  reason,  the  expected  impact  point  is  offset,  the 
problem  of  determining  hit  probabilities  becomes  quite  complex 
Herein  is  developed  a  method  to  approximate  such  a  hit 
probability . 
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I.   INTRODUCTION 

The  probability  of  hitting  a  specified  target  with  a 
given  weapon  system  is  of  particular  interest  to  both  mili- 
tary commanders  and  weapon  systems  analysts.   Determining  hit 
probability  is  not  always  a  simple  matter.   In  particular,  if 
the  impact  points  of  the  weapon  system  are  distributed  accord- 
ing to  some  probability  distribution,  and  the  expected  impact 
point  is  not  at  target  center,  the  problem  of  determining  hit 
probabilities  becomes  quite  complex. 

Over  the  past  thirty  years  this  problem  has  been  studied 
in  depth  as  the  extensive  bibliographies  of  F.  Grubbs  [Ref. 
1] ,  A.  Eckler  [Ref.  2],  and  W.  Guenther  and  P.  Terragno  [Ref. 
3]  illustrate.   Several  approaches  to  solving  the  problem  have 
been  investigated  by  these  and  other  authors  in  an  effort  to 
simplify  the  computations  required  to  estimate  or  actually 
compute  the  hit  probabilities. 

As  an  example  of  one  approach  to  the  problem,  Grubbs  [Ref. 
1]  uses  what  he  calls  a  "single,  straight-forward  and  rather 
simple  technique"  to  approximate  these  hit  probabilities  in 
which  he  employs  an  approximate  central  chi-square  distribu- 
tion with  fractional  degree  of  freedom,  or  a  transformation  to 
approximate  normality.   There  are,  of  course,  other  methods 
available  to  approximate  these  hit  probabilities,  one  of  which 
will  be  discussed  in  the  following  pages. 


The  method  which  will  be  discussed  is  computationally- 
easy,  using  simple  mathematics  in  conjunction  with  a  table  of 
modified  Bessel  functions.   These  tables  are  usually  avail- 
able to  the  analyst  and  enable  him  to  accurately  approximate 
the  desired  hit  probabilities  without  the  use  of  complex  com- 
putations.  It  should  be  mentioned  that  tables  of  offset 
circle  probabilities  have  been  compiled  by  Rand  Corporation 
[Ref.  4].   These  tables  give  hit  probabilities  as  a  function 
of  target  radius  and  offset  distance  of  expected  impact  point, 
but  are  not  as  readily  available  as  the  Bessel  function  tables 
used  in  the  following  procedure.   It  should  also  be  noted 
that  routines  for  Bessel  functions  are  available  for  most  all 
computer  systems. 


II.   NATURE  OF  THE  PROBLEM 


Consider  a  circular  target  with  radius  a.   Place  the 
origin  of  the  coordinate  system  at  the  center  of  the  target 


Let  \  ty.i\\)      be  the  probability  density  function  of  the 
impact  point,  (x,y) ,  of  the  projectile  fired  at  the  target. 
Let  ( f-*>y ,  J^}   be  the  mean  of  the  impact  point  distribution  and 
<J       and  0*  be  the  variance  of  the  impact  point  distribution 

x  *} 

in  the  x  and  y  directions  respectively.  Assume  covariance  is 
zero.  If  the  distribution  of  impact  points  follows  a  bivari- 
ate  normal  distribution,  then 


Assume 


2     9        *> 

0;  =  o;  =  O"   ,  and  R~  h^x  *  Kh    •   Let  ^h  be  the  Proba' 


bility  of  a  hit  on  the  target,  then 


P>T 


y.*y,  -o. 


Changing  to  polar  coordinates  yields 
and 

V1   *       r  -n 


2/hsr 


o      o 


-Jfc7"     *  ,    ~Ztr 


*~  -      r  /Ul  -I 


^^^rc5  Lo 


Rotate  axes  so  that  p.y-R  and  (o.^  -  0    ,    then 


p  .  e^ 

h  " 


*      *  ,.  ** 


,Jr«f(*[j^(t^*)d«]d. 


Looking   at   the    term,      feV-f  ( tS~-  ]d  ©  , 

\co*9d0*  c^TtamO    -(n-iMcos*"  0<nn©(-swO)dO 
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cCn-iMcoC'^sm^ae 

0 


«!  «* 


and  fcM^eae-  (•£)]«*  "ede 
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If   n    is    odd,    then     Cos   Qc\VzQ  .       If   n    is    even,    then 

o 

[cosn0  id  Q  >  (^-0(v\-3WgV--3-3    £^  .      Thus, 

o  A  «... 

UftvccsOWO'^^,  cos  ^e.-^,  ^^Tyr 

-  9-^r  *S  f-\  • -J — -  -»       Also,  for  n  even  and  equal  to  2m, 
*™   -Cj  Vjl)   (yyi\^ 

The  Bessel  function,  denoted  by  JL  ("V)    ,  is  equal  to 

%.  v  ' '  Izj  ,  where  x  is  the  argument  of  the 

Bessel  function  and  p  is  the  order  of  the  Bessel  function. 
The  modified  Bessel  function  of  order  zero  is  denoted  by  ToO0 r 
and  is  equal  to  J0(^\.)  z  f  ^1   s  S   VjJ.  ,  . 

Thus,  fei^l-fccowQldO*  Z*£(!f  'At  =  Z^  T°^> 


*  ah 


so  that  Je*p(~a^sG)dO  =^T0(^) 

It  was  shown  previously  that  the  probability  of  hitting  the 
target,  P   ,  was  equal  to  £^  Lvpts^lp^ l"^^  8l'*,r- 


It  was    also   shown   that   f  e*f>  l^t  cc*  e)  d  Q    -   ZT\  I6  (  ^.  V 
As   a   result,  P  =  £?*      [r£^  ZK  I,  (  %%  )  dr> 


or 


It  is  only  necessary  to  know  the  ratios.  -   and  -  ,  to 

qr     cr 

determine  the  probability  of  hitting  the  target, 

0 

The  Rand  Corporation  has  constructed  tables  giving  these 

hit  probabilities  as  a  function  of  the  two  ratios,  3    and     , 
[Ref.  3]. 

Gilliland,  [Ref.  4],  shows  that 

-&*  "l»  tax1*" 


Luke,     [Ref.    5],    defines    the    function 

jhh)  - 1  -  P^ifc  i&  i'<j  d  t . 

o 

Recall  that    %\    p       yt- 


It  follows  that 


Let  "tr^r"1       dt-  irdir        .   When  t  -  £- „  s  IT ~  <*•  . 


Thus, 

0 

As    a   result, 

h 

Luke,     [Ref.    5],    gives    the    expansion    of     ,J  ( *\t    w\ )         in    a 

series    of    Bessel    functions. 


JK-jVl-e      ^lf\(^  for  v^   >. 

where     ^=     Zlx^'7-     ,  and  Vj  =  ^\  2 . 


Consider 


for    &  *  ^- 
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Then  --l^«1\  - 


for     R.^CL 


where 
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These  latter  two  equations  yield  the  true  hit  probability 
for  a  weapon  system  fired  at  a  circular  target  with  radius  a, 
and  normally  distributed  impact  points  with  expected  impact 
point  at  |Ay  and  JA^  . 

A  good  approximation  to  the  true  hit  probability  can  be 
had  by  considering  only  the  first  few  terms  of  the  infinite 
series  in  the  appropriate  equation.   That  is 


>****\ 


) 


l(MW-e       £ltfwV^        *<*■ 


i.  A*i, 


-•^:)i{%)\m     «>*- 


jh-.i 


where  ir\  is,  hopefully,  relatively  small.   If  v\      is  less 
than  five,  then  the  approximation  would  be  beneficial  to 
analysts  and  would  be  quite  simple  to  use.   If  V\  is 
greater  than  five,  the  calculations  required  to  get  a  reason- 
ably accurate  approximation  to  hit  probability  become 
prohibitive. 
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III.   CONCLUSIONS 
— — . —     ■  ■  ■ 

The  accuracy  of  the  above  approximations  was  checked  on 
an  IBM  360  computer.   The  program  used  the  FORTRAN  language 
and  two  IBM  subroutines,  10  and  INUE,  which  computed  the 
modified  Bessel  functions.   This  program  is  available  through 
the  Operations  Analysis  Department  of  the  Naval  Postgraduate 
School  in  Monterey,  California. 

The  computer  program  was  basically  an  iterative  process 
in  which  different  values  of  the  ratios  ^  and  -  were  set 

o-     cr 

and  the  hit  probabilities  computed  while  the  number  of  terms 

from  the  appropriate  infinite  series  varied  from  one   to  ten. 

The  approximate  hit  probabilities  thus  obtained  were  then 

compared  to  the  true  hit  probabilities  from  the  Rand  tables 

[Ref.  4]. 

The  two  aspects  of  primary  interest  were  (1)  how  accurate 

are  the  approximating  equations  in  the  various  ranges  of  8 

0" 

and  2"  ,  and  (2)  how  many  terms  of  the  infinite  series  in  the 

cr 

approximating  equations  must  be  considered  before  a  reasonable 

degree  of  accuracy  is  obtained.   The  specific  values  of  the 

ratios  3  and  -  ranged  from  .5  to  5.0  in  the  computer  program 
cr     cr 

For  the  equation  for  R<a,  the  ratio  B  was  initialized  at  the 

cr 

value  .5  and  the  ratio  ^  was  allowed  to  vary  between  .5  and 

cr 

5.0.   The  value  for  £  was  then  incremented  to  1.0  and  the 

process  repeated  with  values  of  3-  greater  than  or  equal  to 

cr 

1.0.   This  procedure  was  continued  until  both  B  and  —  were 

cr     cr 
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equal  to  5.0.   For  each  value  of  "  and  -  the  approximate 
hit  probability  was  computed  using  only  the  first  term  of  the 
infinite  series  initially,  and  then  again  for  each  additional 
term  used,  up  to  ten  terms.   In  this  manner  it  was  possible 
to  see  how  rapidly  the  approximating  equations  converged  to 
the  true  hit  probability. 

For  the  equation  for  R>a,  the  same  type  of  procedure  was 
used,  holding  the  ratio  2r  constant  while  the  ratio  6  varied 
through  the  appropriate  values.   As  before,  the  ratio  2-  was 
initialized  at  .5  and  then  incremented  by  .5  each  time  the 
ratio  ?  had  run  its  range  of  values  and  the  hit  probabilities 
computed. 

Representative  values  of  hit  probabilities  obtained 
through  the  approximating  equations  are  contained  in  Appendix 
A.  Figure  1  compares  approximate  hit  probabilities  with  ac- 
tual hit  probabilities  for  R<a,  and  Fig.  2  does  the  same  for 
R>a.  In  both  cases  the  values  shown  were  computed  using  the 
first  five  terms  of  the  infinite  series  from  the  appropriate 
approximating  equation. 

The  convergence  of  the  approximating  equation  values  to 
the  true  hit  probabilities  is  shown  in  the  two  examples  be- 
low.  Example  one  shows  the  convergence  of  the  approximating 
equation  for  R<a  as  the  number  of  terms  from  the  infinite 

series  increases.   Values  of  E  and  2-  were  selected  to  show 

a-  o" 

both  rapid  convergence  and  relatively  slow  convergence.   Ex- 
ample two  shows  convergence  of  the  approximating  equation  for 
R>a.   These  examples  along  with  Appendix  A  point  out  that  the 
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approximation  of  hit  probability  for  R<a  is  optimistic  while 

the  approximation  of  hit  probability  for  R>a  is  conservative. 
Example  1  (R<a) 

*   =  .5,  *  =   1.0 

cr  <r 

Approximation  using: 

one  term  .2089828 

two  terms  .110  870  7 

three  terms  .1047551 

four  terms  .1045011 

five  terms  .1044937       True  probability:   .104491 

six  terms  .1044937 

seven  terms  .10449  37 

eight  terms  .1044937 

nine  terms  .1044937 

ten  terms  .104  49  37 

£   =  4.5,  5*  =   5.0 

cr  o* 

Approximation  using: 

one  term  .9253550 

two  terms  .8596  84  8 

three  terms  .804  4760 

four  terms  .7601165 

five  terms  .7260436        True  probability:   .653154 

six  terms  .7010158 

seven  terms  .6  8342  79 

eight  terms  .6715976 

nine  terms  .6639764 

ten  terms  .6592714 
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Example  2  (R>a) 
a. 

5" 


5    £   = 


1.0 


Approximation  using: 


one  term 
two  terms 
three  terms 
four  terms 
five  terms 
six  terms 
seven  terms 
eight  terms 
nine  terms 
ten  terms 


.0690204 
.0732899 
.0734669 
.0734724 
.0734725 
.0734725 
.0734725 
.0734725 
.0734725 
.0734725 


a.  ft 

cr        <j- 


Approximation  using: 


one  term 
two  terms 
three  terms 
four  terms 
five  terms 
six  terms 
seven  terms 
eight  terms 
nine  terms 
ten  terms 


.0656702 
.1208790 
.1652384 
.1993113 
.2243391 
.2419270 
.2537574 
.2613786 
.2660835 
.2688691 


5.0 


True  probability:  .073473 


True  probability:   .272201 


Accuracy  of  the  approximating  equations  can  be  improved 
in  all  cases  by  considering  more  terms  of  the  appropriate  in- 
finite series.   If  it  is  necessary  to  use  more  than  five 
terms  of  the  series,  however,  the  mathematic  manipulations 
involved  become  impractical  to  do  by  hand.   Accuracy  can  also 
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ft.  CK- 

be  improved  when  -   is  equal  to  -   .   From  Appendix  A  it  is 

0"  ,    <T 

apparent  that  when  =  is  equal  to  ~*  and  both  are  greater 
than  3.0,  accuracy  is  the  poorest.   To  more  accurately  ap- 
proximate the  true  hit.  probability  in  this  case  it  is  suffi- 
cient to  compute  the  arithmetic  mean  of  the  values  obtained 

from  both  equations  (R<a  and  R>a) .   For  example ,  for 

ft  0~  ... 

-       =     -     =5.0,    the    true    probability    is    .460.       The    equation 

cr    cr 

for  R<a  yields  an  approximation  of  .643,  while  the  equation 

for  R>a  yields  an  approximation  of  .325.   The  arithmetic  mean 
of  these  two  values  is  .484,  which  is  more  accurate  than 
either  of  the  two  approximations. 
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IV.   SUMMARY 

If  the  impact  point  of  a  weapon  system  is  distributed  as 
a  bivariate  normal  random  variable,  the  task  of  computing  the 
hit  probability  for  a  circular  target  when  the  expected 
impact  point  is  somewhere  other  than  at  target  center  is  for- 
midable.  The  two  approximating  equations  that  have  been 
developed 

l,  of 


hl.lVl-e^^d^T,^ 


for  R  «£  O. 


te.-o 


for  R>&- 


4**1 

are  relatively  accurate  approximations  to  true  hit  proba- 
bilities using  only  the  first  five  terms  of  the  series. 
Greater  accuracy  can  be  achieved  by  considering  more  terms. 
To  have  an  accurate  approximation  to  such  a  hit  probability 
is  beneficial  to  both  weapon  systems  analysts  and  military 
commanders,  because  of  the  difficulties  involved  in  comput- 
ing actual  hit  probabilities.   The  above  equations  yield  a 
good  approximation  to  hit  probabilities  using  simple  mathe- 
matics and  a  table  of  modified  Bessel  functions,  which  is  not 
difficult  to  obtain. 
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APPENDIX  A 


COMPARISON  OF  HIT  PROBABILITIES 


R    1.0 
-\Tapp   act 


2.0 

app   act 


3.0 

app   act 


4.0 
app   act 


* 


# 


1.0  .267  .267  *     * 

2.0  .731  .731  .410  .397 

3.0  .956  .956  .788  .786   .498  .433     *     * 

4.0  .997  .997  .966  .966   .814  .803   .578  .450 

5.0  .999  .999  .998  .998   .970  .969   .8  35  .813 


Figure  1.   Table  for  R  Less  Than  a 


0. 
ST 


1.0 


R  \  7r    aPP    aCt 


2.0         3.0         4.0 
app   act    app   act    app   act 


1.0  .267  .267 

2.0  .082  .082 

3.0  .011  .011 

4.0  .001  .001 

5.0  .000  .000 


39  3  .39  7 
113  .113 
015  .015 
001  .001 


.399  .433     *     * 
.121  .126   .366  .450 
.016  .017   .120  .133 


5.0 

app   act 


643  .460 


5.0 

app   act 


325  .460 


Figure  2.   Table  for  R  Greater  Than  a 

These  tables  compare  the  approximation  to  the  actual  hit 
probabilities  (app)  and  the  actual  hit  probabilities  (act) 
themselves.   Only  the  first  five  terms  of  the  infinite  series 
described  on  page  11  were  used  for  the  approximations.   Accu- 
racy can  be  improved  in  all  cases  by  considering  more  terms 
of  the  series. 
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